The Stability of Relativistic Neutron Stars in Binary Orbit 
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We analyze the stability of relativistic, quasi-equilibrium binary neutron stars in synchronous 
circular orbit. We explore stability against radial collapse to black holes prior to merger, and against 
orbital plunge. We apply theorems based on turning points along uniformly rotating sequences of 
constant angular momentum and rest mass to locate the onset of secular instabilities. We find that 
inspiraling binary neutron stars are stable against radial collapse to black holes all the way down to 
the innermost stable circular orbit. 
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The two-body problem is one of the outstanding, un- 
solved problems in classical general relativity. However, 
neutron star binary systems are known to exist, even 
within our own galaxy Binary systems are among 
the most promising sources for gravitational wave de- 
tectors now under construction, like LIGO, VIRGO and 
GEO. This has motivated an intense theoretical effort to 
predict the gravitational wave form emitted during the 
inspiral and coalescence of two neutron stars. 

Fully general relativistic treatments of the problem are 
complicated by the nonlinearity of Einstein's equations 
and the need for very large computational resources 
Recently, Wilson and Mathews reported preliminary 
results obtained with a relativistic numerical evolution 
code. Their dynamical calculations suggest that the neu- 
tron stars may collapse to black holes before their orbit 
becomes unstable and the stars plunge. Their results 
are in disagreement with predictions of Newtonian 
post-Newtonian (PN) and perturbation calcula- 

tions 1^, which show that tidal fields stabilize neutron 
stars against radial collapse. 

In a recent paper, we presented quasi-equilibrium, 
polytropic models of fully relativistic, equal mass neu- 
tron star binaries in synchronous circular orbits In 
Newtonian gravity, strict equilibrium exists for two stars 
in circular orbit. In general relativity, because of the 
emission of gravitational waves, binaries cannot be in 
strict equilibrium. However, outside of the innermost 
stable circular orbit (ISCO), the timescale for orbital de- 
cay by radiation is much longer than the orbital period, 
so that the binary can be considered to be in "quasi- 
equilibrium" . This fact allows us to neglect both gravita- 
tional waves and wave-induced deviations from a circular 
orbit to very good approximation. A detailed discussion 
of our approximations and numerical method will be pre- 
sented in a forthcoming paper |p^ . 

In 1^, we focused on the construction of quasi- 
equilibrium models, but did not discuss their stability. 
We found that the maximum allowed mass slightly in- 



creases as the separation of the stars decreases. In this 
Letter we construct sequences of constant angular mo- 
mentum and sequences of constant rest mass, which then 
allows us to apply rigorous theorems based on turning- 
points along the sequences to locate the onset of secular 
instabilities. We find that neutron stars in synchronous 
orbit are stable against radial collapse to black holes until 
they encounter the orbital instability at the ISCO. 

Applying turning-point methods to binaries in corota- 
tion allows us to detect secular instabilities. The separa- 
tion at which simultaneous extrema exist in the mass- 
energy and the angular momentum curves marks the 
point along an evolutionary sequence at which the bi- 
nary becomes unstable in the presence of some dissipa- 
tive mechanism. This instability is distinct from dynami- 
cal instability, which arises independently of any dissipa- 
tion and grows on a dynamical timescale. Gravitational 
radiation can drive the secular instability which is 
expected to occur before (i.e. at larger separation than) 
the onset of a dynamical instability |Q,|2| . It is therefore 
anticipated that it is the secular instabilities that limit 
the range of stable configurations. 

Because of finite numerical resources, we cannot con- 
struct models of very large separation. We therefore ob- 
serve the increase of the maximum mass only for separa- 
tions smaller than about 4 stellar radii. However, even at 
this separation, tidal effects are still fairly small. More- 
over, PN calculations, which are valid for larger separa- 
tions, suggest that the maximum mass should also in- 
crease monotonically with decreasing separation . In 
the following we will therefore assume that the maximum 
allowed mass monotonically increases as the separation 
decreases. 

For a given equation of state, cold equal mass bi- 
nary neutron stars in synchronous circular orbit form a 
two-parameter family, just like single, uniformly rotat- 
ing stars. In our numerical code we adopt the central 
density and relative separation to uniquely specify a par- 
ticular configuration. A Lemma by Friedman, Ipser and 
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Sorkin [|2|, originally derived for single, uniformly rotat- 
ing stars, can therefore be adapted to corotating binaries: 

LEMMA. Consider a two-parameter family of equal 
mass binary stars in synchronous circular orbits based 
on an equation of state of the form P = P(po). Suppose 
that along a continuous sequence of models labeled by a 
parameter X, there is a point Aq at which both Mq = 
dMo/dX and J vanish and where d(f2J + (xMo^/dX ^ 0. 
Then the part of the sequence for which (TlJ -\- fiMo) > 
is unstable for X near Xq . 

Here P is the pressure, po the rest-mass density, J and 
Mq are the angular momentum and rest (baryon) mass 
of one of the two stars, fl the angular velocity and is 
the chemical potential. We will later use M for the total 
mass-energy of one of the two stars (that is half the total 
ADM mass of the spacetime). The proof of the Lemma 
follows directly from Theorem 1 of [|l3|. In practice it is 
more useful to apply the following related Theorem: 

THEOREM 1. Consider a continuous sequence of 
equal mass binary stars in synchronous circular orbits 
based on an equation of state of the form P = P(pq). 
Suppose that the total angular momentum is constant 
along the sequence, and that there is a point Xq where 
Mq = (and where fi > 0, d{fiMQ)/dX ^ 0). Then the 
part of the sequence for which (iMq > is unstable for X 
near Xq. 

The proof follows directly from the Lemma (see 
also [0). 

We can now apply this Theorem to our numerical mod- 
els. In these models the matter obeys a polytropic equa- 
tion of state, P = Kpq'^^^", where K is constant and 
n is the polytropic index. Note that physical dimen- 
sions enter the problem only through K. It is there- 
fore convenient to introduce the dimensionless quanti- 
ties po = K"pQ, Mq = K-"/^Mo, M = K-^/^M and 
J^K-'^J. 

In Figure 1 we plot Mq as a function of the central 
density pc for n = 1. The dotted lines are curves of 
constant relative separation za — ria/r out- Here Tin is 
the inner (coordinate) distance between the two stars, as 
measured between the two closest points along the axis, 
and Tout is the outer distance, as measured between the 
most distant points along the axis. For stars in contact 
Za = 0, and for infinitely separated stars za — ^- The 
lowest dotted line in Figure 1 is for stars at a separation 
Za = 0.3, and the highest curve is for stars in contact. We 
expect that the Oppenheimer-Volkoff curve for infinitely 
separated stars lies below all these lines [Q, such that 
the equilibrium mass supported by a given central density 
monotonically increases with decreasing separation. 

Evolutionary inspiral sequences conserve the rest mass 
Mq and therefore follow horizontal lines in Figure 1. This 
shows, for example, that the central density pQ decreases 
as the two stars approach each other. 

The solid lines in Figure 1 are contours of constant 
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FIG. 1. Contours of constant angular momentum in a 
rest-mass versus central density diagram. The dotted lines 
are sequences of constant relative separation za, ranging from 
0.3 (bottom curve) to (touching stars, top curve). The 
thin solid lines are contours of constant angular momentum, 
ranging from Ji = 0.049 to J2 = 0.0498 in increments of 
0.0001. In the insert we also show contours for J = 0.04955 
and 0.049575. The thick solid lines mark the onset of secular 
instabilities (see text). 

total angular momentum J. The angular momentum of 
infinitely separated binaries is infinite, and decreases for 
approaching binaries. There is a critical angular momen- 
tum J"^*" , for which the contours form a saddle point 
(point A). This critical angular momentum and the crit- 
ical mass Mq'''' at the saddle point take values, which, 
like the maximum mass of nonrotating isolated stars, de- 
pend only on the polytropic index. For n = 1 we find 
J™' ~ 0.05 and M^"' ~ 0.179 

The saddle point divides the turning points of the 
contours into two separate classes: one for contours 
J < J''"', and one for J > J'"''*. We will identify the 
former with the instability at the ISCO, and the later 
with a radial instability of supramassive stars at large 
separations. 

The class for J < J'^'''* exists for all rest masses 
Mq < M"'* and ends at the saddle point A. The physical 
interpretation of these turning points can be understood 
by following an evolutionary sequence for a coalescing 
binary of fixed rest mass, which starts at a large separa- 
tion, e.g., at point F. It evolves along a horizontal line, 
and initially both the angular momentum and the central 
density decrease with the separation. At point E, how- 
ever, the angular momentum goes through a minimum. 
According to the relation 

dM = mj + pdMo, (1) 
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this minimum has to coincide with an extremum in the 
total mass-energy M. This extremum is well understood 
in Newtonian gravity (e.g. Q), PN theory (e.g. 0), as 
well as general relativity and marks the onset of an 
orbital instability at the ISCO. Note that this orbital 
instability is the first instability that these binaries en- 
counter during inspiral. 

For J > J"^^ we find a second class of turning points. 
This class starts at point B in Figure 1 and ends at 
the maximum of the OV curve for isolated, nonrotat- 
ing stars The maximum on this limiting curve 

marks the onset of radial instability for isolated, non- 
rotating stars, and defines their maximum allowed rest 
mass M^^^. By continuity, this class thus marks the on- 
set of radial instability for supramassive sequences, con- 
sisting of stars with rest-masses M^^^ < Mq < Mq"*'. 
Consider, for example, a binary somewhere on the line 
between points C and D. Following the inspiral of these 
stars, they will eventually incounter the ISCO at point 
C. If we could reverse the inspiral and follow an evo- 
lutionary sequence towards larger separations, we know 
that the stars eventually have to become unstable and 
collapse to black holes, since their masses cannot be sup- 
ported when in isolation. This will happen at the turning 
point D. 

Note that so far we have only located turning points 
of contours of constant angular momentum. In order to 
meet all the conditions of Theorem 1 and rigorously es- 
tablish the onset of a secular instability, we also need to 
examine the chemical potential jj,. The chemical poten- 
tial is not normally evaluated by a code that constructs 
equilibrium models (although it could be). Accordingly, 
it is useful to apply a different version of the Theorem 
that does not refer to /i: 

THEOREM 2. Consider a continuous sequence of 
equal mass binary stars in synchronous circular orbits 
based on an equation of state of the form P — P{po). 
Suppose that the rest mass Mq is constant along the se- 
quence, and that there is a point Aq where M — ( and 
where > 0, d{VLM) / d\ ^ Q). Then the part of the 
sequence for which flAI > is unstable for A near Aq. 

The proof follows from the Lemma together with re- 
lation (|l|) (see also [|8|). Note that according to (|l|) the 
extrema in Theorem 1 and 2 have to coincide. 

In Figure 2 we plot contours of constant rest mass Mq 
in an M versus pc diagram. In this plot the contours 
describe secular evolutionary sequences. The contours 
again have two qualitatively different classes of extrema. 
One class are minima of the total energy (eg. point C), 
and the other are maxima (point D). The former can 
again be associated with the onset of an orbital insta- 
bility at the ISCO, while the later mark the onset of 
radial instability of supramassive sequences at large sep- 
arations. Both extrema end at the critical contour M"'', 
above which the contours are monotonically decreasing 
for increasing pc. The value of Mq'''' agrees with the 
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FIG. 2. Contours of constant rest mass in a mass-energy 
versus central density diagram. The dotted lines are se- 
quences of constant relative separation za, ranging from 0.3 
(bottom curve) to (stars in contact, top curve). The thin 
solid lines are contours of constant rest mass, ranging from 
Mo = 0.176 to 0.1794 in increments of 0.0002. The thick solid 
lines mark the onset of secular instabilities (see text). 

height of the saddle point A and point B in Figure 1. All 
extrema in Figures 1 and 2 agree to within a few percent, 
which is a measure of the accuracy of our numerical code. 

In order to verify the assumptions in Theorem 2, we 
have to make sure that the extrema of the constant rest 
mass contours do not coincide with an extremum of the 
angular velocity. In Figure 3 we show the two contours 
connecting points E and F, and C and D, together with 
corresponding plots of f2. 

The minima at E and C can by analyzed very easily. 
Using central density as the parameter A, all the assump- 
tions, M = 0, > and d{(lM)/d\ ^ are satisfied, 
and we conclude that the configurations to the left of the 
minima, for which f]M > 0, are unstable. These minima 
mark the onset of a secular instability close to the ISCO. 

The maximum at D is a little less obvious from our 
numerical data, since it occurs very close to a minimum 
in Cl. However, the offset is measureable and larger than 
our numerical error, and we consistently find this offset. 
In addition, we know that supramassive stars must be- 
come unstable when moved out from finite separation to 
infinite separation. We therefore conclude that the ex- 
trema in M and are indeed distinct, so that we can 
apply Theorem 2. This establishes the onset of a sec- 
ular radial instability for the maxima on the thick line 
extending downward from point B in Figures 1 and 2. 

We conclude that all binary configurations in the area 
under the thick solid lines and dashes in Figures 1 and 2 
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at smaller separation than secular instabilities, together 
with the absence of a secular instability against radial 
collapse, up to the ISCO, suggests that a dynamical in- 
stability against radial collapse may also be absent in 
normal binary neutron stars. 

This work was supported by NSF Grant AST 96-18524 
and NASA Grant NAG 5-3420 at lUinois, and NSF Grant 
PHY 94-08378 at Cornell. 
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FIG. 3. Contours of constant rest mass. The top panels 
show two contours from Figure 2, and the bottom panels show 
the corresponding angular velocities. 



are stable. This area can be subdivided into a normal and 
a supramassive region. All sequences with rest masses be- 
low the maximum allowed rest mass M^^^ of an isolated, 
nonrotating star are normal sequences. These sequences 
start at finite separation (eg. F) and end at the ISCO 
(eg. E). In addition to these normal sequences, there are 
supramassive sequences with rest masses between M^'^^ 
and M"'*. Evolutionary curves for these sequences con- 
nect a radial instability at large separation (eg. D) with 
an orbital instability at the ISCO (eg. C). We conclude 
that all normal binary neutron stars in synchronous or- 
bit, with rest masses below the maximum for isolated, 
nonrotating configurations, are secularly stable against 
radial collapse to black holes down to the ISCO. At the 
ISCO they encounter a secular orbital instability. 

Supramassive stars are unstable to collapse at large 
separation, but, if formed at finite separation, remain 
stable as they inspiral down to the ISCO. Supramassive 
stars could arise in principle if a rapidly rotating, supra- 
massive core were to undergo bifurcation at finite radius. 
Alternatively, they could arise if two spinning neutron 
stars become bound in a binary, or if two neutron stars 
accrete matter while in a binary. 

The secular instability evolves on the same timescale 
as the inspiral, since both are driven by gravitational ra- 
diation P,pT[. Orbital plunge will occur at the onset of 
dynamical instability, which defines the true ISCO. More- 
over, realistic binaries are more likely to be closer to irro- 
tational (circulation = 0) than corotational jl9j . For such 
systems, however, the location of the dynamical instabil- 
ity is quite close to the onset of the secular instability, 
both in Newtonian theory and in PN theory . The 
onset of the dynamical instability is at slightly smaller 
separation. 

The fact that dynamical instabilities typically arise 
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